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The fluctuations of 7 ray families simulated with CORSIKA in the energy region 3x10^^^- 
lO^'^eV on the basis of standard Lns collider physics exhibits alignments of secondaries in the 
stratosphere and at ground level. The remarkable event registrated on the Concorde doesn't 
fit well however those cases ; The possible hints of new mechanisms, especially the valence 
diquark breaking, are considered. Observing that the extrapolation of the original cosmic ray 
primary spectrum derived from the size spectrum measured in the Akeno classical EAS array 
coincides with the spectrum measured recently by the Hires Stereo experiment, we point out 
a possible overestimation of the primary energy in inclined showers of the surface arrays like 
AGASA. 



1 Remarkable events with coplanar emission 

The attention on coplanar emission was motivated by the events recorded in Pamir X ray 
chamber'^', in Kambala, in Kascade, but also in balloon experiments and in the low stratosphere 
with high resolution X ray emulsion chambers in the Concorde^^ the geometrical criteria used 
to select an alignment treat directly the coordinates of the individual 7's, either the linear 
coefficiently 
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(1) 



or the parameter [1] defined as: 
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(2) 



where (p'^j is the angle between the straight Hnes joining the i^^ and j^^ particles to the A;*^ one 



(0 < (p'j^j < vr). Alignments have been easily obtained, independantly on the models (as far as 
they have similar transverse momenta distribution), and practically on the primary mass and 
energyEl. 



2 The Concorde event near 10 eV 



Tracing back the aligned events produced with CORSIKA at Concorde altitude, we have ob- 
served that an unbalanced pt received in the first collision on an energetic vr'^ could produce 
an alignment when the cascade starts 10km above the chamber. Among the 211 7's of the 
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Table 1: Calculated fractions (%) of the aligned events with at least 4 7-rays (7+e*) above 10 TeV for with 
jrj > 0.94 5(first row) and A4 > 0.8 (second row) for different high energy hadronic interaction models. 

JF2aF2 Concorde event, we show on a lego plot the energy deposited in a plane perpendicular 
to the axis by 34 7's (about one half of the total energy deposited, i.e. 1600 TeV). Inside the 




Figure 1: Left: Lego plot of the energy deposited by 34 7's in the alignment (energy in TeV on vertical axis, 
X, y in mm). Right: An example of the histogram of invariant mass for cluster A, with the mass of the tt'^, the 

maximum gives an interaction distance near 100m 

alignment, the four most energetic 7's with respective energies, 300 TeV, 105 TeV, 75 TeV and 
53 TeV are lying on a perfect straight lineal with r = 0.9993. The invariant mass histograms 
for the different clusters (Fig^ suggest an interaction level at about 100m above the Concorde, 
in contradiction with the simulation on the basis of standard physics requiring about 10km. 
Such circumstance suggests a violent phenomena separating the valence quarks of the proton 



projectile, especially the valence diquark, which cannot be recombined with a quark of the sea, 
suppressing the leading baryon at the end of the collision. 

One hypothesis is that the original rotating relativistic string between the valence quark 
and the valence diquark becomes a more complex system with a secondary string (centered on 
the barycenter of the diquark) between the two quarks partners of the diquark. The maximal 
tension of the strings occurs when the quarks are at the largest distance from each other, i.e. 
when the 3 partons are on a common diameter which would be the axis of the fragmentation. 
The shorter mean free path of the diquark in the target nuclei could help the phenomena. As 
a consequence of the suppression of the leading particle, the maximum depth Tmax is expected 
to level off during one energy decade above the kneeS and some typical behavior would have to 
be observed in EAS ( in muon electron dependance, in the age parameter versus size, different 
absorption length, enhanced steepeness of the most energetic 7's and hadrons spectra). Prom 
other part, there have been few simulations with primary nuclei, and the recent proposition of 
a predominant a component suggests to examine more carefully the asymmetries in light nuclei 
collision, i.e. effects like the giant dipole resonance. 



3 Lateral distributions in giant EAS 

Some functions are used in large surface arrays without reference to the total size or to the age 
parameter, giving just an interpolation between the detectors to evaluate the densities required 
for the estimators at 600 m or 1000 m from axis. The couple {N(.,s) is especially useful in the 
case of hybrid events^l at the level of the registration; it can be derived from the fluorescence 
measurements to start a minimisation on the densities recorded with the surface array and give 
a better determination of the axis position with the hypergeometric function^ 

f{x) = g{s) x'-^ix + 1)^-^(1 + dx)-" (3) 

with a = 1.92, b = 3.8, c = 7.71, d = 0.00342 (for distances r and densities A, r = xvq and 
A = T-/(r),ro = 36.8m, N being the total size. The calculation of the HG serial is replaced by 

the approximation valid up to 6=40° : g{s) = — 0.19 + 0.969 s — 0.468 s^. This relation 
works also for vme's with a = 1.94,6 = 3.92, c = 2.87, d = 0.00562, tq = 39.2m, g(s) being 
replaced by gj,me{s) = 1 - 0.789 s + 0.133^^ and finally Ayme= 4>i{Ne/rl)f{x) with = 0.47 
(agreement with experimental data in ref [6]). 



4 Primary spectra from classical and giant arrays 

The good extrapolation of the spectrum obtained in Akeno with the spectrum from HIRES 
Stereo is shown in Figure 2. In Akeno^^, the densities were determined with a modest detector 
spacing (30m or 100m) and a specific lateral distribution, containing the age parameter was 
employed to localize the core and obtain the size N . The size N is converted directly to the 
primary energy with a relation in agreement with CORSIKA within 2%. The 20 km^ array 
(Array 20) with 19 detectors, separated by about 1 km from each other, uses the distribution: 
p{r) = N Ce x~°'{l + x)~(''~")(l + 2M))~'^^ ^^'^ normalisation constant). This analytic descrip- 
tion with a fixed value a = 1.2, without reference to the age parameter is used to determine the 
axis position and to interpolate the value of the density at 600m. 

In contrast to the size conversion in Akeno, the scintillator response in terms of density Seoo 
is here converted to the primary energy following: E2o{eV) = 2.010^^ x (^eoo)^'" (S(r) is related 
to the electron and muon densities) 

In place of the size spectrum, the 5*600 differential spectrum in Array 20 is obtained taking 
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Figure 2: DifTerential primary spectrum for Array l(full square), Array 20(open square), AGASA(full triangle) 
and HIRES Stereo experiments(open triangle). Fits to Akeno(full line) and Array 20(daslied line) are froml^ for 
the clarity of the graph, the error bars are not plotted for AGASA data. 

an attenuation length Agoo in parallel to Ag in Array 1 following: 



A constant value (also for AGASA) A6oo= 500(7 ~ cm ^ was employed. The most recent values 
reported by AGASA^^ are more close from the values of Akeno than the values of Array 20 
(figure 12) ; the intensities of AGASA remain however larger than for Array 1 in the overlapping 
energy region and exhibit a general excess by 30% when compared to Hires Stereo data^^. Prom 
our simulation data, we have derived the values of the attenuation length Agoo for different 
zenith angles: for small inclinations Q < 30° the values of the attenuation length concerning 
proton primaries are quite more important than the average value A6oo= 500g' — cm~^ used in 
AGASA. When the primary energy is increasing, the depth of the maximum becomes more and 
more close of the arrays in altitude, such as AUGER or AGASA : the conversion of inclined 
densities to ^goolO) according to equation |1] becomes poorly appropriate as the cascade is far 
from a stable absorption phase, especially for protons primaries. In the depth interval of about 5 
radiation units following the maximum, the absorption process is described by the age parameter 
increasing in parallel from 1.0 up to 1.2, the lateral distribution around 600m from the axis 
becoming flatter. The increase of this flattening of the density distribution turns to a systematic 
overestimation (via the vertical density from relation(4), the shower recorded may be classified in 
bins of larger energy). Above 3.5 lO^^eV, a clear divergence in the discrepancies between AGASA 
and Hires Stereo appears rising from 150% above 300% at 6.10^^ eV. This may come again from 
the lateral distribution becoming flatter more rapidly than the reduction of the total size : the net 
result is that the densities (at 600 m) are 5 — 10% larger in the bin G = 20° — 30° than the vertical 
density when the atmospheric depth separating the array and the shower maximum becomes 
lower than 3 cascade units. Some systematic errors could also enter in the axis localisation. 

5 Conclusions 



S'6oo(0) = 5'6oo(0) X exp{ 




(4) 



A large proportion of the alignments may be explained by fluctuations, however, the alignments 
observed in the stratosphere indicate the necessity of a more carefuU analysis and the collection 
of new events, in the low stratosphere. 5000 Hours could be available for a scientific payload 



during the certification of the Airbus A380. The flights carried at an altitude of 13.1 km 
[nOg/cm"'^) with 10 emulsion chambers, similar to those used in Concorde, would multiply by 
100 the statistics of remarkable 7 ray families. This remains the most simple approach to the 
behavior of the valence quarks at energies close of the LHC energy range. The present approach 
points out a better consistency between the spectra obtained by classical size measurements and 
Hires Stereo measurements, favourable to the GZK prediction. The spectrum measured by the 
array KASCADE-Grande will be useful to improve the calibration of giant surface arrays. 
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